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We examine the N = 2 Wess-Zumino model defined on the d = 2 Euchdean 
lattice in connection with a restoration of the Leibniz rule in the limit a ^ in 
perturbatively finite theory. We are interested in the difference between the Wilson 
and Ginsparg- Wilson fermions and in the effects of extra interactions introduced 
' by an analysis of Nicolai mapping. As for the Wilson fermion, it induces a linear 

^ ■ divergence to individual tadpole diagrams in the limit a ^ 0, which is absent 

in the Ginsparg- Wilson fermion. This divergence suggests that a careful choice 
' of lattice regularization is required in a reliable numerical simulation. As for the 

■ effects of the extra couplings introduced by an analysis of Nicolai mapping, the extra 

, couplings do not completely remedy the supersymmetry breaking in correlation 

I functions induced by the failure of the Leibniz rule in perturbation theory, though 

those couplings ensure the vanishing of vacuum energy arising from disconnected 
diagrams. Supersymmetry in correlation functions is recovered in the limit a — > 
QhI with or without those extra couplings. In the context of lattice theory, it may be 

properly said that supersymmetry does not improve ultraviolet properties but rather 
it is well accommodated in theories with good ultraviolet properties. 



X. 

S ■ 1 Introduction 

It is important to define the supersymmetric Wess-Zumino model on 
the lattice in such a way that the non-renormalization theorem [)2| Q and 
consequently the absence of quadratic divergence is preserved. Since 
the absence of quadratic divergence arises from a subtle cancellation of 
bosonic and fermionic contributions, we have to ensure the precise (not 
approximate) bose-fermi cancellation. This task is not easy if one recalls 
that the Leibniz rule is generally broken on the latticej^]. See Refs.0- 
TJ] for the analyses of related issues. 



Recently, it was suggested [lT5j] that a perturbatively finite theory, if 



latticized, could preserve supersymmetry to all orders in perturbation 
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theory in the sense that the supersymmetry breaking terms induced by 
the failure of the Leibniz rule become irrelevant in the limit a ^ 0. It was 
demonstrated in Ref . [|T5[] that this is in fact realized if one first renders the 
4-dimensional Wess-Zumino model finite by applying the higher deriva- 
tive regularization. A non-perturbative confirmation of this proposal has 
not been given yet, but we believe that a perturbative confirmation of the 
absence of quadratic divergence is a prerequisite for the non-perturbative 
analysis^. In Ref.[jrB|], the Ginsparg- Wilson fermion [[T71] [[T^j [[T9|] was uti- 
lized, which has a nice chiral property but at the same time introduces 
certain subtle aspects to the analysis [Til . 



A 2-dimensional reduction of the Wess-Zumino model, which exhibits 
N = 2 supersymmetry, is finite perturbatively with qualifications to be 
specified later, and thus it provides a good testing ground of the sugges- 
tion made in Ref. [jT5|, though the crucial issue of quadratic divergence 
cannot be studied in this model. In the present paper, we examine the 
N = 2 Wess-Zumino model on the d = 2 Euclidean lattice perturbatively 
and clarify several basic issues involved which may become relevant in an 
actual numerical simulation. Firstly, we examine the use of the Wilson 
fermion instead of the Ginsparg- Wilson fermion since the Wilson fermion 
is much easier to handle numerically. However, the Wilson fermion in- 
duces a strong chiral symmetry breaking and thus it is important to see 
if this introduces any new aspect into the problem. The suggestion in 
Ref. [[T5|] is based on making all the Feynman diagrams finite, namely, the 
cancellation of divergences among Feynman diagrams is not sufficient in 
general. If one applies this criterion to the present problem, we encounter 
one-loop level divergences in some of the individual Feynman diagrams for 
correlation functions though those divergences cancel among bosonic and 
fermionic contributions. (In disconnected vacuum diagrams, two-loop 
diagrams contain divergences.) In particular, the Wilson fermion intro- 
duces a linear divergence to individual tadpole diagrams in d = 2, which 
is absent in the Ginsparg- Wilson fermion. The presence of linear diver- 
gences suggests that the lattice regularization is not arbitrary but need 
to be "well-behaved" one, which ensures the precise cancellation of these 
linear divergences among diagrams for a 7^ 0, such as in a formulation 



^It should be noted that the lattice in this context is introduced not to control the divergences but to 
make numerical and other non-perturbative analyses possible. 
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with precise lattice supersymmetry in the free part of the LagrangianH. 
If all the Feynraan diagrams should be absolutely convergent, the latti- 
cization would enjoy more freedom to recover supersymmetry in the limit 

The second issue analyzed is the role played by extra couplings in- 
troduced by an analysis of Nicolai mapping [j5j. The Nicolai mapping 
in the present context suggests an appearance of exta interactions [fT]] [pTT 



[|T^ which vanish if the Leibniz rule is satisfied on the lattice. Also these 
extra terms spoil the naive hypercubic symmetry on the lattice. Because 
of these novel features of the extra interactions, one may hope that these 
extra terms might remedy the failure of the Leibniz rule appearing in 
the remaining conventional interaction terms. We analyze this issue in 
the framework of perturbation theory. Our result shows that these extra 
terms do not completely remedy supersymmetry in correlation functions, 
which is broken by the failure of the Leibniz rule, at least in weak cou- 
pling perturbation theory, though those terms ensure the vanishing of 
vacuum energy arising from disconnected diagrams. As far as the corre- 
lation functions are concerned, supersymmetry is recovered in the limit 
a ^ with or without those extra couplings. 



2 N=2 Wess-Zumino model and Nicolai mapping 

We start with a Lagrangian defined in terms of the Wilson fermion on 
the d = 2 Euclidean lattice 

jC = ^(L'(i) + L'(2))^-m^^-2^V^(P+(/)P+ + P_(/)*P_)^ 

+ ra(2)c/)+(ra(2)0r 

+ #2(Vf + ^Vf)0 + ^(02(vf + iVf)c/)r (2.1) 

where i/j is a two-dimensional Dirac spinor. Since we are interested in 
the d = 2 model as a toy model for 4-dimensional theory, we choose the 
superpotential to be a specific form 

W'{(J)) = m0 + g(j)\ (2.2) 

^The actual numerical analysis, for example, is simplest in the simplest form of the Lagrangian, but 
due care is required to ensure supersymmetry in the limit a 0. 
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Here we defined 

= ^—{ilj{x + all)^ijj{x-afi)-2i/;{x)). (2.3) 

We note the important property Y.x f{x){V'^g){x) = — Ex(V^/)(x)p(x). 
Our Euclidean 7 matrix convention is 

75 = 75, 

^± = ^(1±75). (2.4) 

When we have the operator in the bosonic sector, we adopt the 

convention to discard the 2x2 unit matrix. The terms 

Ck^n = - + F^F (2.5) 

stand for the kinetic (Kahler) terms. The last two terms in (2.1) are the 
extra terms introduced by an argument of Nicolai mapping[[7|, while other 
terms are the naive lattice translation of the continuum theory except for 
the Wilson term and its super-partners. Namely, the terms 

jCw = ipD^2)i/^ + ^(2)0 + {FD^^2)4>T (2.6) 

stand for a naive supersymmetrization of the Wilson term, which induce 
a hard breaking of continuum chiral symmetry. 

The last two terms in (2.1) vanish if the Leibniz rule for Vj^ should 
hold on the lattice, namely, if (tP'[x){y^^(t))[x) = \{y^^(l)'^){x). This fact 
suggests that the extra terms might remedy the supersymmetry breaking 
induced by the failure of the Leibniz rule in the remaining interaction 
terms. These extra terms also break the (hyper-)cubic symmetry on the 
lattice. 

The elimination of the auxiliary fields F and in the starting La- 
grangian gives 

C = ^(L'(i) + L'(2))^-V5(P+H^" + P_(H^")> 
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ay 

+ W(Vf + zVf)(/)+ (W(Vf + iVf)0)* 
with W = m(J) + 5f0^ by noting 

E V^0(a;) = - E V^0(a;)c/)(a;) = 0. 

X X 

This Lagrangian agrees with the one introduced in Refs. 
analysis of Nicolai mappingi. 

The Nicolai mapping here is defined by 

6 = (Vf + L'(2))A-[/-Vf5, 



6 = (-Vf + L'(2))5-l^-VfA 



s 



(2.7) 

(2.8) 
n|] |T2| by an 



(2.9) 



with 



(t){x] 
U = 



V2 



(A + iS), 



1 



y = -^xw'-iwr). 

If one uses the specific representation of 7 matrices 



7 



1 

-1 



7 



-1 
-1 



75 



i 
-i 



(2.10) 



(2.11) 



the Jacobian for the transformation from (^1, ^2) to (A,B) precisely agrees 
with the determinant of the fermion operator in (2.7). The bosonic part 
of the Lagrangian is then written as 

E Ao.on(a:) = -T.l[^l{x)+&{x)] (2.12) 



and the partition function is given by 



Z 



Vi/jVi/jVAVB exp[E C{x)] 



= /%2^6exp{-E|[e?(2^)+^l(^)]} (2.13) 



^Recently, the Nicolai mapping was extended to the d = 2 Wess-Zumino model defined in terms of 
Ginsparg- Wilson operators, which makes chiral symmetry manifest . 



if one imposes universal (periodic) boundary conditions both on fermionic 
and bosonic variables. The vanishing of the vacuum energy is thus en- 
sured by the Nicolai mapping even for a 7^ 0. The partition function is 
reduced to 



if the continuum limit a ^ is well-defined. The presence of the Nico- 
lai mapping would then ensure the degeneracy of bosonic and fermionic 
spectra of H in the continuum limit. 

3 Supersymmetry transformation 

One may define a lattice supersymmetry transformation parametrized by 
a constant Dirac-type Grassmann parameter e by 



The supersymmetry transformation parametrized by a constant Dirac- 
type Grassmann parameter e, which is treated to be independent of e, is 
given by 



Tr(-l)^exp[-/?^] 



(2.14) 



6^ = -e[P_0 + P+0*]i?(i) - e[P_F* + P+F] 

Si/j = 

6^ = eP+tp, 6(J)* = eP-iJj 

5F = eP_D^i)ip, 6F* = eP+D^i^ip. 



(3.1) 



6i; = -P>(i)[P-0 + P+01e - [P-F* + P+F]e 

S(j) = ijjP+e, 6(1)* = 'ijjP-e 

SF = ij)D^i)P_e, 6F* = ipD^i^P+e. 



(3.2) 



Here we note that 




(3.3) 
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Under the transformation (3.1), it can be confirmed that the kinetic term 

J6Ck^n = J{-e[P^(l) + P+(l)^]D^i)-e[P^F^^P+F]}D^,)ilj 
+ / (l)*d,d,eP+ilj + / (t)d,d,eP^ilj 

+ / F^eP^D^i)ilj + / FeP+D^i^i/j = (3.4) 

is in fact invariant. The mass term 

^mass = —mtjjtjj — mF(/) — mF*(/)* (3-5) 

and the Wilson term 

jCw = V5i?(2)^ + ra(2)0 + (ra(2)0)* (3.6) 

are also confirmed to be invariant under the above supersymmetry trans- 
formation by using the relation D(i)D(2) = F){2)D{i)- 
The variation of the (conventional) interaction terms 

Cnt = -2gtjj{P+(l)P+ + P^(t>^P^)ilj - g[F4>'' + [Fcf'f] (3.7) 

is given by 

jdC.nt = g J e[-2P4D^,)cl))cl)i/j - 2P4D^i)cl)*)(l)*iJj 

+P4L'(i)02)^ + P+(L'(i)(0*)2)V;] (3.8) 

where we used the relation Ex / (3:) (i^(i)fi') (a;) = —J2x{D{i)f){x)g{x). This 
variation of the interaction terms would vanish if the difference operator 
Z)(i) should satisfy the Leibniz rule 2{D(^i^(J)){x)(J){x) = {D(^i^(J)^){x). (The 
terms quadratic in i/j vanish by themselves.) 

As for the U{1) x Uji{l) charges (and holomorphicity) analogous to 
those in the 4-dimensional theory [p], we may assign 

0=(1,1), 
^=(1,-1), 

P+^ = (1,0), P_V^ = (-1,0), 
m= (-2,0), 

^ = (-3,-1) (3.9) 
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for the terms appearing in the conventional formulation. Even for m = 
= 0, the Wilson term Cw in the Lagrangian violates the U (1) symmetry. 
As for the extra terms introduced by an argument of Nicolai mapping 

jCextra = g(J)\Vf + iVl)cJ) ^ igcj)\Vf + iViW (3.10) 

which break hypercubic symmetry, they are not invariant under the above 
supersymmetry transformation. Since the lattice supersymmetry trans- 
formation we defined above preserves transformation properties under the 
hypercubic symmetry, the supersymmetry variation of these extra terms 
do not mix with the variation of other terms. Those extra terms also 
break Ur{1) symmetry. We re-iterate that these extra terms vanish if the 
Leibniz rule should hold on the lattice. 



4 Feynman rules for perturbative calculations 

It is interesting to examine to what extent the extra terms introduced by 
an argument of Nicolai mapping [j7|[lnj][|T^, namely the last two terms in 



C (2.1), are essential to maintain supersymmetry in perturbation theory. 
Although the extra terms, which break hypercubic symmetry, do not 
mix with other terms under the lattice supsersymmetry transformation in 
(3.1) as we already explained, the lattice supersymmetry transformation 
is not unique and thus we cannot a priori exclude the possible cancellation 
of supersymmetry breaking effects among the interaction terms@. Our 
assumption is that perturbative calculations are universal at least for 
a small coupling constant and independent of the specific definitions of 
lattice supersymmetry transformation. One starts with the free part of 
the Lagrangian 

Co = tjj{D^i^ + D^2))'ip-(l)*Dl^D^,)(J)^F*F-rnijtj;-rn[F(l)+{F(l))^ 

+ FL'(2)0 + (ra(2)c/))^ (4.1) 

When we have the operator Z)|^<)i?(i) in the bosonic sector, we adopt the 
convention to discard the 2x2 unit matrix. 



"'in fact, an exact Ward identity which is regarded as a part of supersymmetry is known to exist in 
the Lagrangian defined by the Nicolai mapping [pi] . We shall analyze this identity later. 
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The propagators are given by 



L'|i)L'(i) + (-m + L'(2))2' 
(FF*) = {-)—, , 

= (^Vt) = (-) W^^f^f^^ 

^(1)^(1) + (-"^ + ^2))^ 

W = 7^ = - ^ 4-2 

and other propagators vanish. Here we have Z)|-^^ = — -D(i). In the mo- 
mentum representation, we have 



(FF*) = (-) 



(sin a/c^)^ 



(F0) = (F>*) = 



(sinaA;^)2 + (aM)2' 
a(aA1) 



where 



- ai^i^sinak^^ a{aM) . . 

= (smafc,)^ + (aM)^ 

7W(a/c^) = - cosa/c^) + m. (4.4) 

The interaction terms for perturbative calculations are given by 

Cint = -2g^{P+(jyP+ + P-(j)'P-)'il^-g[F(j)'' + {Fct>y] 

+ gcf'iVf + iWDc^ + gic^^Vf-riVlW. (4.5) 

If one sets one of in the extra interaction terms (i.e., in the last two 
terms in Cmt) to be a constant 

(j){x) = 00 = constant (4-6) 

then the extra terms go away by noting Ea; 4>{x)'^^i4>{^) = 0. This means 
that the effects of the extra terms introduced by the Nicolai mapping do 

9 



not appear in the one-loop diagrams if one sets the momenta of external 
lines at 0. In other words, only those diagrams where the momenta of 
external lines cannot be set to be zero are affected by those extra terms. 



5 Lower order diagrams 

We now examine several lower order diagrams for correlation functions in 
perturbation theory^. The theory in d = 2 becomes more convergent in 
higher order diagrams, and one can confirm that the possible supersym- 
metry breaking effects in higher order diagrams become irrelevant in the 
limit a ^ 0, provided that one-loop sub-diagrams are properly treated. 
The one-loop diagrams are thus crucial in the analysis of supersymmetry. 
As for the disconnected vacuum diagrams, they shall be later analyzed 
separately. 



5.1 Tadpole diagrams 

The one-loop tadpole diagrams for the scalar consist of two diagrams; 
the first one is a scalar loop and the second is a fermion loop. The scalar 
loop contribution is given by 



-2g(l){F(p) = -2g(t) / tt^-T- 1 \2 ^ ( ka\ 



2 



J-n {27rY^J {sink^Y + {aM{k^)Y ^ ^ 
where we re-scaled the integration variable in the last expression as 

ak^ k^ (5.2) 

and defined 

aA1(/c^) = ~ cos /c^) + am. (5.3) 

In the limit a ^ this diagram diverges as ~ 1/a, namely, linearly 
divergent. This divergence is worse than the divergence in continuum 



perturbative analysis of the 4-dimensional Wess-Zumino model with the Wilson fermion was per- 
formed in 1^. In the present d = 2 model, the perturbation is based on g/m ^ 1. 
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theory (and also in the lattice theory with the Ginsparg- Wilson fermioni) , 
which is logarithmic. 

The fermion loop contribution is given by 

which is precisely cancelled by the scalar contribution (5.1) for a finite 
a. However, each diagram is linearly divergent due to the strong chiral 
symmetry breaking by the Wilson term. In a numerical simulation, one 
would need to choose the free part of the Lagrangian to be lattice su- 
persymmetric, as in the present formulation, so that the cancellation of 
linear divergence is exact^. Alternatively, one may introduce an auxiliary 
regularization such as higher derivative regularization to make each dia- 
gram convergent and thus less sensitive to the parameter a in the limit 
a ^ 0, as we discussed in 4-dimensional theory [|T^ . 

5.2 Induced 0^-coupling 

We have contributions from a scalar loop and a fermion loop. The scalar 
loop contribution is given by 

^i2gf<p{FcP){F<p)<P 

_ 9 2/2 d'^k a{aM) a{aM) 

9 ^ i-TT/a (27r)2sin2(a/c^ + ap^) + {aMf {sinak^f + {aMf 

(fk {aM{k^ + ap^)) 



2^ H-Pfi)^iPf^) 



^ (27r)2 sin {ak^ + ap^) + {aM{k^ + ap^))' 



{s\nk,f {aM{k,)f ^ ■ ^ 

which approaches a constant for a ^ 0. This behavior is consistent with 
the continuum behavior, but the difference is that all the momentum 



^The power counting with the Ginsparg- Wilson fermion is identical to that of continuum theory p^. 
^If all the Feynman diagrams should be absolutely convergent, one would enjoy more freedom in 
choosing lattice regularization. I thank H. Kawai and T. Onogi for a helpful comment on this point. 
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regions, not only the infrared region, contribute to the integral. This is 
an effect of the chiral symmetry breaking by the Wilson term. 
The fermion loop contribution gives 



2! 

- -2a^M-v Uiv ) r {aM{k, + ap,)) 

{s\nk,f + {aM{k,)f ^^-""^ 

which is precisely cancelled by the scalar contribution (5.5) even for a 
finite lattice spacing a. 

Each Feynman diagram which gives induced couplings higher powers in 
(j) such as (j)^ is reduced to the continuum result in the limit a — ^ 0. Those 
couplings in any case cancel among the scalar and fermion contributions 
even for finite a. The non-renormalization of the superpotential in this 
sense is thus maintained in the one-loop level, and in higher-loop levels 
in the limit a ^ 0. 

5.3 Self-energy corrections 

The simplest self-energy correction is that to the auxiliary fields F and 
F*. The one-loop correction is given by 



= 2/FF^(00*>(00^) 



2 



9 _ f-TT/a (Pk c? a- 

= 2a FF / 

^ i-TT/a (27r)2 sin^iak^ + ap^) + {aMf (sin ak^f + {aMf 

2 7-./ \T7*r r 1 

= 2g F[p,)F [p,] ^2ii)Hm^k,^{aM{k,)f 

X — (5 7) 

sm^{k^^ap^)^{aM{k^^ap^)f 

This integral vanishes for a — ^ if one keeps the integration domain 
outside 

\k^\ < 6 for all fi (5.8) 
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for arbitrarily small but finite 5 and for fixed p^. The integration inside 
the above domain gives a finite continuum result if one notes 

(sin k^f + {aM{k^)f ^k^^ {^k^ + amf ^kl + amkl + ("^)^ (^•^) 

inside the above domain. A re-scaling of k^ back to the original momen- 
tum variables k^ ak^ gives the continuum result in the limit a ^ 0. 

The fermion self-energy correction is given by 

i(2^)2^(P+0P+ + P_0*P_)V^V^(P+0P+ + P_0^P_)^ 

^ (2^)2[v;P+(V^V^)P_^(#*)+V5P_(V^V^)P+^(#*)] (5.10) 



7r/a d k ai^^ sm{ap^ ^ ak/j) 



a 



{2g) i/^iPf,) J_^^^ {27ryKm^ (ap^^ak^) + (aM)2 (sin a/c^)2 + (a>()2^'^^^^^ 

- ^9^2 7/ . d^k aiY sm{ap^ ^ k^) 1 

- [^9) ^Wy_,(27r)2^sin2(ap^ + /c^) + (aA^)2(sin/c^)2 + (a>()2J'^^^^^- 

This integral vanishes if one sets j)^ = 0, which means that the fermion 
mass receives no quantum correction when renormalized at vanishing 
momentum, despite the chiral symmetry breaking by the Wilson termi. 
This integral also vanishes for the domain outside (5.8) in the limit a ^ 0, 
and the integral is reduced to the continuum result for the domain inside 
(5.8) in the limit a ^ for fixed p^. 

To analyze the wave function renormalization, we consider the case 
with an infinitesimal p^^ but the lattice spacing a kept fixedly. Namely, 

Ip^I <C m, 1/a. (5.11) 

We thus expand 

sm{ap^ -\- k^) ~ sin k^ + ap^ cos /c^, 
aM {k^ + ap^) ~ aM (k^) + ap^, sin k^, 
sin^iap^ + k^) + {aMf ~ sin2(A;^) + {aM{k^)f 

-^ap^ sin 2k^ + 2aM.{k^)ap^ sin /c^. (5.12) 

^This vanishing mass correction arises from the differences of the Feynman rules in the present model 
and QCD. Also, all the higher loop corrections are reduced to the (supersymmetric) continuum results 
in the limit a — > in the present model. The Wilson term does not always imply the mass shift. 

^It should be noted that we assume a small coupling g/m <C 1 and infinitesimally small external 
momentum but otherwise make no assumtion about the lattice spacing a. 
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The integral in this expression is given by 

2 /-TT (fk S/i i^Pf^ cos kf^ 

Y^ii ^7'^ sin k^,[py sin 2kv + 2{aM{k))pT^ sin ky] 
[s\n\k^) + {aM{k)YY 



} 



2 (fk E/z i7>M cos /c^ 



-7r(27r)2 ^[sin^(/c^) + (a7W(fc))2]2 



1 . . 5 1 

+ - E sin ^/^EP.^[g.^2^^^)^^^^^^))2]2} 



^ 1 2 r ^^fe r i7^PM cos fe^ , . . 

2 i-TT (27r)2 ^ [sin2(/c^,) + {aM{k)yY^' ^ ^ 

By noting the symmetry under ki ^ k2, we thus have the wave function 
renormahzation for the fermion 

which disagrees with the finite renormahzation factor for the fields F and 
F* at = in (5.7) for a finite a; 



(27r)2 [(sin/cJ2 + (a;t^(/,^))2]2 (27r)2 [sin2(/c^) + (a>t(/c))2]2 

TT (fk a^lT,„{l — cosk„) 



d^k o? j-TT (Pk a2|EzyCos/Cj, 

2 
1 

- (27r)2 [(sinA;^)2 + {aM{k^)y]^ ^ ^ ^ 

for a 7^ 0, though this difference vanishes in the limit a — ^ 0. This shows 
that the finite wave function renormahzation factor breaks supersymme- 
try for a 0. 

We next examine the self-energy corrections to the scalar field (j). The 
contribution from a scalar loop diagram in the conventional interaction 
terms gives 

l^g^lFcj^^ + (Fct^^ym^ ^ (Fct^^y] 

-^g\{Fct?fFct?) 

= Ag^\F^F){<^^c^)c^ (5.16) 
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J TT 



X 



(sinfcj2 



(27r)^ 



(sin A;^)2 + {aMf sin2(A;^ + aip^) + (a>l(/c^ + a^?^))^ 
The one-loop fermion contribution is given by 



-{2gf4>''Tr[{P+^i^P^){P^ijijP^)]c^ 



=-(2.)m: 



2 ./-Tr (27r)2 

^ r Z7^ sin(fc^ + api^) i-f^' sin k^ 

^ ^^sin2(/c^ + ap/,) + {aM{k^ + ap^))2 (sin/c^)2 + {aMf^^ 



£k 



(27r) 



2 



r sinfa + aj)^)sinfc^ 1 

^sin2(/c^ + ap^) + (aAi(/c^ + apj)2 (sin k^y + (aA^)2J- ^ " ^ 

The sum of these two contributions gives rise to 



(27r)2 

(sin(fc^ + gp^) - sin k^) sin fc^ 1 

''W(/c^ + apJ + (a7W(/c^ + ap^))2(sinfc^)2+(a7W)2J ^ 

which vanishes for = 0. This means that the mass correction to the 
scalar particles exactly vanishes in the one-loop level. However, each 
term logarithmically diverges in the limit a ^ 0, which suggests that 
the choice of the free-part of the Lagrangian should be at least invariant 
under the lattice supersymmetry transformation to ensure the divergence 
cancellation, such as in the present formulation. This integral vanishes for 
a ^ for the domain outside (5.8) and for fixed p^j^. For the the domain 
inside (5.8) and for fixed p^, the integral is reduced to the continuum 
result in the limit a ^ 0. 
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For an infinitesimal p^, we have 
4g (f) [Pi^mp^,) / TTT^I- 



(27r)2^ (siny2 + (a7W)2 
^ _d_ 1 _ |E^[(Qj3^)^sin^A:j 



= 4^/^0*(p^)0(p^) ^{Y,ap^ cos k^ sink ^ 

1 ^ 1 |EM[(Q^PM)'sin^fe/^] 

2 [sin2(A;^) + {aM{k^)YY [{smk^f + (aA^)2]2-^ 

/■TT d^k 

= -2g^(f)*ip^)(f)ip^) ^^^2 Y^IpI cos 2A;^ + sin2 fcj 



X 



[sin2(/c^) + (aA^(/c^))2]2 
= -2g^*ip,)plcl>ip,) ^a'^ E[cos 2fc^ + ^(1 - cos 2k,)] 

[8in\k,) + iaM{k,))Y ^^-^^^ 
This deviates from the renormahzation of F and F* for finite a, 

i-TT (fk 

J-^ (27r)2[sin2/c^ + (aM(A;J)2]2 

TT ^2/. Q,2i ^^[^Qg2A:j, + i(l - cos2/c^)] 
-vr (27r)2 [sin2 k,^{aM{k)Y]^ 
^ (fk a^\Y.u{l - C0S2K) 

(27r)2 [sin2 k, + (aA^(/c^))2]2 ^ 

though this difference vanishes in the limit a — >^ 0. 

5.4 Self-energy corrections induced by extra couplings 

Finally, we examine the effects of the exta couplings in (4.5) introduced 
by an argument of Nicolai mapping on the self-energy of scalar particles. 
This is given by 

^[gct>\Vl + iVf)0 + g{ct>\Vl + zVf 
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- /[(0'(Vf + ^Vf)c/))1[02(Vf + iVf )0] 
2/((Vf + ^Vf )0)^(0*0)(0»(Vf + iV?)0 

+V((vf + ^vf )0)*(0*0)(r(vf + ^vf 

+V^*^((Vf + zVf)0)'^c/))(0»(Vf + zVf)0 
+4gYmf + iVf)0)»(0*(Vf + iVf)0)<^ 
+V(/>*(((Vf + iVf)0r(Vf + iVf)0) (0*0)0. (5.21) 

The first term in (5.21) gives 

2/((vf + ivimp,)nvf + ivi)<p(p,) ^ 

1 



sin2(/c^ + ap^) + (a7W(/c^ + a]?^))2 sm^{k,,) + (a7W(/c^))2 

X I — (5 22) 

sin^ {kf, + ap^) + {aM{kf, + ap^))^ sin2(fc^) + {aM{kf,)y'^ ' 

This gives for an infinitesimal 

The second term gives 

, 9,,.sinapi sina»2N ,/ nn^,/ ^ d'^k 

a -i sm{ki + api) + sin(A;2 + ap2) . . 

sin^ (kf,) + {aM{kfj,)y sin^ik^ + ap,,) + (a>l(/c/, + a]?/,))^ 

which gives for an infinitesimal 

d?k a 

— 7Vjtii pns Z:", -I- /TDr, pns Z:".^ 



•^-TT (27r)^sm (/c^) + (aAl(/c^))"' 

za]?i cos ki + a]?2 cos A;2 



sm\k^) + (a>l(A;^))2 
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d 1 

+ (-z sin ki + sin /C2) E aPi^^i — • 2/, n , / a^^/ ^^2 ] 

-TT (27r)^ sin^(A;^) + {aM{k^)y 

1 . —iapi cos A;i + ap2 cos k2-. 
""2^ sin^(/c^) + (aA^(/c^))2 ^ 

= 9^{-m - P2)(t>{Pny(l>{Pti) I !t \. • 2/7 X , ^/ ^^9 

•^-T^ (27r)2sin^(/c^) + (a7W(A;^))2 

(-ia]?i + aj)2) Ez/Cos/Cj, 
sin^(/c^) + (aM(fc^))2 ^ 



The third term gives the complex conjugate of the second term, which 
agrees with the second term itself. Thus we have altogether 

r> 2 i / \~k 2 1 / \ J2v COS ki, /r r>/^\ 

-2s 4>ip,) P,^M (2^)2[,i„.(fc^) + (^^(fc^))2]2 - (5-26) 

The fourth and the fifth terms give precisely the negative of the contribu- 
tions of the conventional interactions to the self-energy of scalar particles, 
(5.16) and (5.17), which we have already evaluated. 

If we collect all the terms arising from the extra couplings together, 
we obtain 

^ 2i/ 2i/ \ d^^k a^[l - Ez/cos/cz. + iEi.cos2/cj, + iEi.(l - cos2A;yl]^_, 
2fy{p,) P^0(P.)/_, (2^)1 [sin^(M + {aMikJw ^''^ 

which vanishes in the limit a ^ 0, as it should be since the conventional 
interaction already ensures supersymmetry in the limit a ^ 0. These 
terms do not help the wave function renormalization factor of (p agree 
with that of either F or ijj. The breaking of supersymmetry in the wave 
function renormalization factors persist for a 7^ even if one includes the 
effects of the extra couplings induced by an analysis of Nicolai mapping. 
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6 Low-energy effective action with one-loop correc- 
tions 



The low-energy effective action which includes the one-loop quantum cor- 
rections is written in a momentum representation as 

,2 



771/ 

Ceff = (1 + fk/j - mil;^ - (1 + z^)(f)*pl(f) - -— — (f)*(f) + .. 

J- ~r Zp 



— 777, _ 

= (1 + z^)[ipi - V^V^] 

J- ~r 



(6.1) 



after the elimination of the auxiliary fields F and F*. Here we defined the 
finite wave function renormalization factors (see (5.7), (5.14) and (5.19)) 

^ J-^ (27r)2[(sin/c^)2+(aA^(A:J)2]2' 
_ 2 r d'^k I a'^^Eiy cos ki, 

_ 2 r a^[^Ez.cos2A:^ + |Ezy(l - cos2A;^) 

" ^ J--J^^ [sin^(y + (aM(/c))2]2 J" ^^"''^ 

Supersymmetry suggests the uniform wave function renormalization 

1 + 2;^ = 1 + 2;^ (6.3) 
and the degeneracy of the mass parameter 



2 2 

m m 



(6.4) 



{l + z^y {1 + zf){1 + z^ 
or in the accuracy of one-loop correction 

2z^ = {zp^ z^). (6.5) 

If one includes the contributions from the extra couplings, these condi- 
tions are replaced by 

Zip Z(i) -\- Z^xtrai 

2z^ = {Zp + Zcl, + Zextra) (6.6) 
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with (see (5.27)) 

_ 2 r a^[— 1 + Ejycos/cj, — |Ei/Cos2/ci, — i Ezy(l — cos2A;j.)] 

= 2z^ - zp- z^. (6.7) 

It is interesting that the degeneracy of the mass parameter, namely, the 
second condition in (6.6) is satisfied even for finite a in the presence of 
the extra couphngs. However, the uniform wave function renormahzation 
condition 

= zf (6.8) 

is still broken since z^ < zp for finite a. The supersymmetry is thus 
broken for a 7^ even with the extra couplings induced by the Nicolai 
mapping. 

In the continuum limit a ^ 0, we have 

Z^p = Zp = Z(j), 

Zextra = (6.9) 

and the supersymmetry is recovered, with or without the extra couplings. 
This conclusion is valid up to any finite order in perturbation theory. 

7 Check of Ward identity 

The Nicolai mapping suggests that the Lagrangian is written as 

£ = ^''i^)^^/iy) - 1 uu^)r (7.1) 

where^ 

{A^} = {AB), (7.2) 
and thus we have the relation 

- {r(^)r(v)} = I T>i^^eM-T.\j:iU^)n (7.3) 



^°We identify the spinor index of the Dirac fermion with the flavor index of the scalar particle, an 
apparently Lorentz non-invariant operation. 



20 



which is equal to 



as can be confirmed by expanding Aa{x) formally in powers of ^^{z). 
These relations give rise to the identity [jn|] [jT2] 



{r{xW{y)) + {Aa{x)ip{y)) = 0. (7.5) 

We check this identity for a small momentum region. The fermion 
propagator with one-loop quantum corrections is given by 



(^^) = = 1 + toli + Oipl) (7.6) 



in the low-energy limit | ;^/m| <C 1 but with fixed a. 
We next note 

U = mA^^{A'-B'), 



V = mB + ^AB. (7.7) 



We evaluate 

75 



{U^)Aiy)) = iVt^D^,^){Aix)Aiy))-m{Aix)Aiy)) 

^ {{A' - B'){x)A{y)) - Vl{B{x)A{y)) (7.8) 



V2 

with the interaction terms 



Cnt = -^{{m-D^,))A){A'-B')-^{{m-D^))B){AB) 



:(V? A - V^B){A' - B') + -^{-VfB - V^A){AB) 



'^^ tJ;{A^iB-fr,)tJ;. (7.9) 



V2 
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We first have 



(V? + D^,)){Aix)A{y)) - m{A{x)A{y)) 
—ipi — m 

" (1 + + ^ 

= - - '11^^ + 0(pI). (7.10) 



m m 

We next evaluate 



^ {[A'-B')(x)A(y)). (7.11) 



V2 

The contributions from the conventional interaction terms give in mo- 
mentum representation 

^ (27r)2 (sin/cJ2 + (aM(A;^))2 (sin/c^ + ap^f + (aM(/c^ + ap^)f^ 



X 



(sinap^)2 + {aM{ap^)y 

(7.12) 

which gives for small p^ 

^. (7.13) 

m 

The contributions from the extra couplings give in momentum represen- 
tation 

2 r F 1 

^ (27r)2 (sin/c^)2 + (aMik^))^ (sin/c^ + a]?^)2 + {aMik^ + apj)^-^ 
ai sin api 



(sinap^)2 + (aA1(ap^))2 
2//"^ c/^fc a i(sin /ci -|- api) 

~ ^ (27r)2 (sin + {aM{k^)y (sin A;^ + ap^f + + apj)^ ^ 



a' 



XT-zrvrTTrTJ7zr^2 (7-14) 



(sinap^)2 + {aM{ap,j)) 
give for small 

,2 



(zi^-z^)^. (7.15) 
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This term vanishes for a ^ 0. 

These calculations show that the Ward identity [|TTj] []T^ 

{tPtP)u + {Ux)A{y))=0 (7.16) 

is precisely satisfied up to the order 0{pj^), 



\2h 



Ux)Aiy)) = iVf + D^,)){A{x)Aiy))-m{Aix)A{y)) 



V2 



^ {{A' - B'){x)Aiy)) -Vl{B{x)A{y)) 



(I + zf) ipi{l ^ zf) zf , Ap\ 



V{zf- z^)^ + 0{p 



m m'^ m ^ 

= -^-'^^±^ + 0(pl) (7.17) 

even for 7^ at a 7^ 0, if one recalls (|^)ii = Pi- But in any case the 
correction terms induced by the extra interactions vanish ^zf — z^ ^ 
in the limit a ^ 0. Consequently, the Ward identity in the limit a ^ is 
not modified by the extra interactions introduced by the Nicolai mapping. 
This is consistent with the numerical findings of the behavior of various 
exact and not-exact Ward identities, which appear to be equally valid 
numerically in the limit a ^ [O . 



8 Disconnected vacuum diagrams 

Contributions to the vacuum energy cancel exactly for the free part of the 
Lagrangian due to the precise lattice supersymmetry. The lowest non- 
trivial contributions from interaction terms arise in the two-loop level. 
The (conventional) two-loop scalar contribution is given by 

^g\F4>' + {Fct?y][F4>'' + {Fct>y] 
4 g\[F4>'yFct>') 

= 2g^{F^F){cl,^ct>){ct>^ct>) (8.1) 

(i^P fn £k 1 



(27r)2 J-n (27r)2 (sinp^)2 + iaM{p^)y 
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{smk^f 1 

The two-loop contribution which contains a fermion loop is given by 

^ (2^/)VP-0*P-V^V^P+0P+V^ 

^ -(2^)2(0»rr[(P+^^P_)(P_^^P+)] 

" 2 J-n (2^ Itt (27r)2 (sin p^^ + (a7W(p^))2 

.. y^r i7^sin(/c^+j3,,) i-f^" sin k^ 

^sm\k^ + j)^) + {aMik^ + ]?^))2 (sin k^y + (aA^(]?^))2^^ 



4^ 



0^2 J) rn crk 



(27r)2 ^-TT (27r)2 {sinp^^f + (aA^(]?^))^ 



r sin(/c^+j3^)sinfe^ 1 

^sin2(A;^ + {aM{k,^p,)Y {smk.y + (aA1(/c^))2J- ^ " ^ 

In the limit a ^ 0, these two integrals contain logarithmically divergent 
infrared sigularities in the re-scaled variables p^ and k^. These diver- 
gences, which agree with the divergences in continuum theory, precisely 
cancel each other. However, the remaining finite parts of these two inte- 
grals do not quite cancel each other even in the limit a ^ and thus lead 
to the non- vanishing vacuum energy. This is a result of supersymmetry 
breaking by the failure of the Leibniz rule. This complication arises since 
the vacuum diagrams are not finite even in d = 2 (in fact contain loga- 
rithmic overlap-divergence) and all the momentum regions contribute to 
these vacuum diagrams. 

A way to remove these finite contributions in the limit a ^ (without 
relying on the extra couplings) is to apply a higher derivative regularization 
which amounts to the replacement of all the terms in the free part of the 
Lagrangian (4.1) as 

£0 = ^(P'(i)+P'(2))i?^-0"P'{i)P'(l)i?0 + P*i?P 

- mtjjRtjj - m[FR(J) + (FRcJ))*] + PP>(2)i?(/) + {FD^2)R(J)T (8.3) 
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where R is the higher derivative regulator 

^ m ^^■^> 

with a new fixed mass scale M. This regularization preserves the lattice 
supersymmetry (3.1) and (3.2) in the free part of the Lagrangian. By this 
way, all the non-vanishing contributions are limited to the momentum 
regions pj^ < M^, and the vacuum diagrams comletely cancel in the limit 

It is interesting to see how the extra couplings help to remove the 
vacuum energy even for a 7^ 0. This is given by 

i[^02(Vf + iVf )(/> + ^(0^(Vf + iViW] 

X [#2(vf + ,vf )0 + ^(02 (Vf + iVlW] 

= V(((vf + ivf)0)»(0*(vf + ivf)<^)((/>» 

+2/(((Vf + iViWiVf + iVf)0)(0»((^» 

^ (2^ (27r)2 (sinp^)2 + (aM(pJ)2 

r sin(fc^ +p^)sin fc^ 1 

^^sin2(/c^ + ]?J + {aM{k^ + p^)y {sink^y + (a>l(A;^))2^ 

(fk 1 



+2/ r 

J — TT 



(27r)2 (27r)2 (sinp^)^ + {aM{p,)f 

i (8 5) 

(sin k^y + {aM{k^)y sm^k/, + p^) + {aM{kf, + p^))2 



We thus confirm that the vacuum energies (8.1), (8.2) and (8.5) put 
together completely cancel for finite a, and this is a nice aspect of the 
analysis based on the Nicolai mapping. 

9 Discussion and conclusion 

We have examined the N = 2 Wess-Zumino model on d = 2 Euclidean 
lattice in connection with a restoration of the Leibniz rule in the limit 
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a ^ 0. In particular, we examined the Wilson fermion instead of the 
Ginsparg- Wilson fermion. We also examined the effects of extra couplings 
introduced by an analysis of Nicolai mapping. 

As for the Wilson fermion, it introduces linear and logarithmic di- 
vergences in some of individual Feynman diagrams, though those diver- 
gences precisely cancel among Feynman diagrams for correlation func- 
tions in the formulation which ensures supersymmetry for the free-part 
of the Lagrangian. In the general analysis of the Leibniz rule in Ref . [[TB| 



each Feynman diagram was made finite to ensure the Leibniz rule in the 
limit a ^ 0. In such a case, the lattice regularization would enjoy more 
freedom since it is introduced just to allow the numerical and other non- 
perturbative analyses, and the lattice artifact is safely removed in the 
limit a ^ 0. Each Feynman diagram in the N = 2 Wess-Zumino model 
in d = 2 which was analyzed here, however, is not finite in general in 
particular with the Wilson fermion, and the precise cancellation of these 
divergences for finite a is important. 

As for the effects of the extra couplings introduced by an analysis 
of Nicolai mapping, which breaks hypercubic symmetry, these couplings 
do not completely remedy the breaking of supersymmetry induced by 
the failure of the Leibniz rule, though those extra couplings ensure the 
vanishing vacuum energy. We also illustrated how the Ward identity [jn|] 
[[T2|] is satisfied even if the uniform wave function renormalization, which 
is required by supersymmetry, is not satisfied for finite a. 

For a minimal latticization of the Wess-Zumino model in d = 2 which 
ensures lattice supersymmetry for the free part of the Lagrangian but 
without the extra couplings 

jC = tjj{D(^,)^D(^2))'^-mtj;tlj-2gtj;{P+(l)P+ + P^(l)*P^)tlj 

+ FL'(2)0+(FL>(2)0r, (9.1) 

we have illustrated that all the supersymmetry breaking terms in correla- 
tion functions induced by the failure of the Leibniz rule are irrelevant in 
the sense that those terms all vanish in the limit a ^ 0. This is consistent 
with the general analysis of perturbatively finite theory on the lattice [[TS 
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The lattice operation implies 



{V{fg)){x) = {Vf){x)g{x) + f{x){Vg){x) + a{V f){x){Vg){x) (9.2) 

if one defines {Vf){x) = {f{x^a) — f{x))/a, and thus the breaking of the 
Leibniz rule is of order 0(1) if the momentum carried by field variables 
is of order 0(l/a). To the extent that the derivative of field variables 
is required in supersymmetry to balance the dimensionality of bosonic 
and fermionic variables, the Leibniz rule is indispensable for the validity 
of supersymmetry. It is well known that supersymmetry improves ultra- 
violet properties of field theory. In the context of lattice theory, one may 
rather reverse the argument and one may even argue that the finite theory 
is required to accommodate supersymmetry in a consistent manner since 
the conventional definition of derivative 

dfjx) ^ ^.^ f{x + g) - f{x) 
dx a^o a 

which satisfies the Leibniz rule presumes that the momentum carried by 
the field variable f{x) is negligibly small compared to 1/a. This is realized 
in lattice theory only if the theory is finite at least in perturbative sense. 

In conclusion, our analysis is consistent with the past analyses of the 
d = 2 Wess-Zumino model and we believe that our analysis gives an expla- 
nation why these past non-perturbative numerical analyses worked [|nj] [jr2[ 



in particular, both of the Ward identity which is expected to be exact 
on the lattice and those Ward identities which are expected to be broken 
by the lattice artifacts0[|r2j]. A numerical calculation of the mass cor- 



rection also appears to be consistent with the (continuum) perturbative 
estimate, as was noted in [|T^. All the supersymmetry breaking effects 



in correlation functions induced by the failure of the Leibniz rule become 
irrelevant in the limit a ^ for a finite theory. The existence of the 
Nicolai mapping in the d = 2 Wess-Zumino model is a nice property of a 
specific formulation of the specific model, such as ensuring the vanishing 
vacuum energy, but it is not a necessary condition for a consistent defi- 
nition of supersymmetric models on the lattice in the limit a ^ 0. The 
finiteness is a more universal condition which ensures supersymmetry in 
the limit a ^ 0. 



^^I thank S. Catterall for a helpful communication related to this issue. 
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Finally, the analyses of other aspects of supersymmetry on the lattice, 
which were not discussed in the present paper, are found in Refs.[|2D|]-[|231 . 



Note added: 

M. Golterman and D.N. Petcher[[2^|] analyzed related issues in the 
context of N = 1 Wess-Zumino model in d = 2. I thank M. Golterman for 
calling the above work to my attention. S. Elitzur and A, Schwimmer [|25| 
discussed a related problem in a Hamiltonian formalism. I thank A. 
Schwimmer for calling their work to my attention. 
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